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Abstract. We review the extension of the factorization formalism of perturbative QCD to coherent soft 
rescattering associated with hard scattering in high energy nuclear collisions. We emphasize the ability 
to quantify high order corrections and the predictive power of factorization approach in terms of uni- 
versal nonperturbative matrix elements. Although coherent rescattering effects are power suppressed by 
hard scales of the scattering, they are enhanced by the nuclear size and could play an important role in 
understanding the novel nuclear dependence observed in high energy nuclear collisions. 

PACS. 12.38.Bx - 12. 39. St - 24.85. +p 



1 Introduction 

Rescattering in hadron-nucleus and heavy-ion collisions 
provides an excellent tool to diagnose properties of nuclear 
medium and could play an important role in understand- 
ing novel nuclear dependence recently observed at rela- 
tivistic heavy ion collider (RHIC) and in planning future 
experiments at the Large Hadron Collider (LHC). Many 
approaches in studying the rescattering effects have been 
proposed and used for calculating nuclear dependence in 
high energy nuclear collisions [1,2,3 ,4 ,5 EE). 

In this talk, we focus on a treatment of coherent soft 
rescattering associated with hard probes |2EH- Our work 
is based on perturbative QCD (pQCD) factorization ap- 
proach, which is different from the works of Baier et al. 
(BDMPS) [1] and Zakharov [H], and the reaction opera- 
tor approach of Gyulassy et al. 0. The BDMPS analysis 
does not require the presence of a hard scattering, but de- 
scribes the coherent results of many soft scatterings. Its 
primary subject is induced energy loss. Our analysis re- 
quires a hard scale, and begins with the pQCD treatment 
of hard-scattering with emphases on momentum transfer, 
caused by coherent initial- and final-state soft scatterings 
2,3. Our work attempts to stay as close as possible to the 
pQCD factorization formalism, in which we may readily 
quantify high order corrections in powers of strong cou- 
pling constant a s , as well as corrections that decrease with 
extra powers of momentum transfer |5l5tll)lll| . 

In the following we consider only initial- and final-state 
interaction that gives leading power in medium length 
[A 1 / 3 ) and in a s at each scattering. We first identify the 
coherence length in nuclear collisions and the source of the 
leading medium size enhancement. We then apply pQCD 
factorization approach to calculate the leading nuclear de- 
pendence in several physical observables. We show that if 
we neglect soft rescattering off quark fields |12U13| . the 



Fig. 1. Sketch for the scattering of a parton (or a lepton) in a 
large nucleus with a momentum transfer . 



leading medium effects induced by multiple soft rescat- 
tering depend on only one well-defined nonperturbative 
matrix element, (F +a F a + ), defined below. We extract its 
value from different physical measurements and discuss its 
universality. A brief summary is given at the end. 



2 Coherent multiple scattering and leading 
nuclear A-dependence 

A hard probe corresponds to a scattering process with a 
large momentum exchange whose invariant mass Q = 



q*\ » ^Iqcd, as sketched in Fig. ^ It can probe a 
space-time dimension much smaller than a nucleon at rest, 
1/Q -C 2i? ~ fm, with nucleon radius R. But, the same 
probe might cover a whole Lorentz contracted large nu- 
cleus, if 1/Q > 2R(m/p) with averaged nucleon momen- 
tum p and mass to, or equivalently, x -C x c = \/2mR ~ 
0.1 with x being an active parton momentum fraction in 
the scattering, xp ~ Q. The critical value x c corresponds 
to the nucleon size. If the active x is much smaller than 
x c , a hard probe could cover several nucleons in a Lorentz 
contracted large nucleus and interact with partons from 
different nucleons coherently. 

Inclusive deeply inelastic scattering (DIS) on a nucleus 
offers an ideal example of coherent multiple scattering 
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Fig. 2. Coherent multiple scattering of the struck parton 

in deeply inelastic scattering (a) and in t-channel in hadron- Fig. 3. (a) Diagram with poles that give rise to an yl 1/3 -type 
nucleus collisions (b). enhancement to DIS cross section off a large nucleus; and (b) 

corresponding twist-4 quark-gluon correlation function. 



and power corrections The strength of the scatter- 
ing is defined by the virtual photon momentum, q^. Let 
q^ = —xbP^ + Q 2 /{2xbP ■ n)n^ with defined along 
a lightcone direction opposite to that of p^ . While the 
Q 2 = —q 2 sets up the hard collision scale, the scattered 
quark probes the nuclear matter via multiple soft final- 
state interactions. When Bjorken xb = Q 2 /2p • q <§; x c , 
the multiple scatterings at a given impact parameter are 
coherent over entire size of the Lorentz contracted nu- 
cleus. This is best seen in the Breit frame, as shown in 
Fig.[2Ia), where the incoming quark reverses its direction 
after interacting with the virtual photon and collides with 
the "remnants" of the nucleus at the same impact pa- 
rameter. The same coherent multiple scatterings can take 
place in hadron-nucleus collisions along the direction of 
momentum exchange of the scattering |14j. as shown in 
Fig.Hb). 

To identify the leading medium length enhanced nu- 
clear effects, let's consider multiple scattering contribu- 
tion to inclusive DIS off a large nucleus, as sketched in 
Fig- Eta). For a spin-averaged inclusive DIS cross section, 
there is only one large momentum scale, Q, "invariant 
mass" of the exchanged virtual photon, and the factoriza- 
tion is expected to be valid for all power corrections as a 
consequence of operator product expansion (OPE) [T5]: 



da DIS = d °2 c 

Q 2 



1 + 



z^a 2 



Tr 



i/A 



dot 

Q 4 



l + C ( 1 - 4 >a s+C ( 2 'V + 
l + c ( 1 ' 6) a s W 2 ' f V + 



■i /. 1 



T, 



i/A 



(1) 



where <g) represents convolutions in fractional momenta 
carried by partons and T„ represents a parton correla- 
tion function or a matrix element of a twist n operator. 
In Eq. the o % n and c"' n J are perturbatively calculable 
short-distance partonic cross sections or coefficient func- 
tions, which are independent of the target size. There- 
fore, we need to find the nuclear size (A 1 / 3 -type) enhance- 
ment induced by multiple rescattering from the matrix 
elements, if there is any. 

The leading twist parton distributions, T 2 s, represent 
probability densities to find a single parton in a target and 
can have some nuclear dependence via the input distribu- 
tions to their DGLAP evolution equations. The A x / 3 -type 



target-size enhancement to DIS cross section can only ap- 
pear in the terms beyond the first row in Eq. • For defi- 
niteness, let's consider a leading power suppressed contri- 
bution to DIS cross section, as shown in Fig.[3Ja), which 
can be factorized into the form according to Eq. Q 2 , 



da A = X! / dXl dX2 dX3 T (ii')/PA ( Xl > X2 ' x ^ 
fan J 



(«') 
x do 



( 4 ) I \ 

l i i, ) (xi,x 2 ,x 3 ) . 



(2) 



The matrix element , as sketched in Fig. EJb) , is typ- 
ically of the form [5], 

T(ii>)/ PA (x 1 ,x 2 ,x 3 ) 
dy^dy 2 dy^ 



i ip + (x 1 y 1 +x 2 y 2 +x 3 y 3 ) 

(2tt) 3 " 

x (pA\BH0)Blfo)B v fo)Bite)\pA) , (3) 

where p = pa/A and Bi is the field corresponding to a 
parton of type i = q,q, G. The structure of the target is 
manifest only in the matrix clement T in Eq. J2Jl. Each 
pair of fields in the matrix element Eq. © represents a 
parton that participates in the hard scattering. The y~ 
integrals cover the distance between the positions of these 
particles along the path of the outgoing scattered quark. 
In Eq. |J3J| , integrals over the distances y ~ generally cannot 
grow with the size of the target because of oscillations from 
the exponential factors e %p XiV * . 

Since the kinematics of a single-scale hard collision is 
only sensitive to the total momentum from the target, two 
of the three momentum fractions: x\,x 2 , and X3 cannot 
be fixed by the hard collisions. Therefore, there is always 
a subset of Feynman diagrams, like one in Fig. Ofa) with 
poles labeled by the crosses "x", whose contribution to 
the partonic parts, in Eq. © is dominated by regions 
where two of the three momentum fractions vanish. The 
convolution over dx\ dx 2 dx^ in Eq. J2J) is simplified to an 
integration over only one momentum fraction, 

f 3 

J II dXt T (jj')/pA(i x i}) tfj3>)({ x i}p) 

dxT q (x,A)o[^{xp), (4) 
where the partonic part i s 

finite and perturbative, 
with the superscript (D) indicating the contribution from 
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(a) 



(b) 




Fig. 4. Tree diagrams give to the leading medium-size en- 
hanced contributions to DIS cross sections. 



double scattering. The above matrix element, T q (x, A), as 
illustrated in Fig. EJb) , has the form 

T q (x,A) = J d £e^ J d A^0 { y- 2 - yr)% 3-) 
x l(pA\q(Oh + F a+ (y^)F+ a (y^) q (y^)\p A ) (5) 

where \pa) is the relevant nuclear state. The variable x 
here is the fractional momentum associated with the hard 
parton from the target that initiates the process. Similar 
gluon-gluon correlation function T g (x,A) is also impor- 
tant In this type of the twist-4 parton-parton cor- 
relation functions, two integrals over the y~ and y 2 can 
grow with the nuclear radius. However, if we require local 
color confinement, the difference between the light-cone 
coordinates of the two field strengths should be limited 
to the nucleon size and only one of the two y~ integrals 
can be extended to the size of nuclear target. The twist-4 
parton-parton correlation functions are then proportional 
to the size of the target, that is, enhanced by A 1 / 3 . 

It is important to emphasize that using a pole in the 
complex Xi (longitudinal momentum) space to do the in- 
tegral does not assume on-shell propagation for the scat- 
tered quark. Indeed, the Xi integrals are not pinched be- 
tween coalescing singularities at such points, and the same 
results could be derived by performing the Xi integrals 
without going through the Xi — points 0. 



3 Coherent power corrections in DIS 

The inclusive DIS cross section on a nucleus provides an 
unique probe for effects of coherent multiple scatterings by 
varying the value of Bjorken xb- Under the approximation 
of one-photon exchange, the unpolarized inclusive DIS 
cross section probes two structure functions: Ft(xb,Q 2 ) 
and Fl(xb,Q 2 ), corresponding to the transverse and lon- 
gitudinal polarization states of the virtual photon, respec- 
tively |17II18| . The structure functions at the lowest order 
in a s are given by ^7] 

F^ T \x B ,Q 2 ) = jE e !^'^) + %) , (6) 



F^ T '(x B ,Q*)=0(a s ) , 



(7) 



where (LT) indicates the leading twist contribution - the 
first row in Eq. p]l. runs over the (anti)quark flavors, 



e q is their fractional charge, and q(x,Q 2 ) is the leading 
twist quark distribution: 

Q 2 ) = J % e ixp+v ~ <Pl?(0) y q(y-)\ P ) (8) 

in the lightconc A + = n^A^ = gauge for hadron momen- 
tum p» = p+n^ with = [1, 0, 0_J and n» = [0, 1, 0J. 

The Feynman diagrams in Fig. 0] give the leading tree 
level contributions to the lepton-nucleus DIS cross section. 
The cut-line represents the final state [3]. For transversely 
polarized photons Fig. Ufa) gives the leading twist par- 
tonic contribution, 

1 

2 



da^ = ^e 2 q 5(x 



xb) , 



(9) 



from which F^ LT ' ) in Eq. |JBJ is derived after convoluting 
with the leading twist quark distribution in Eq. |j8J). Dia- 
grams with two gluons in Fig. 01b) gives the first power 
correction to transverse partonic cross section |3lll8| , 



da, 



(i) 



X -e 2 
2 q 



2N r 



-,2 r 



2tt 2 F 2 (0) 



xb 



dx 



5{x 



with the two-gluon field operator 
dy^dyl 



xb) 
(10) 



F\0) 



(27T) 



F +a (y2)F a + (y^)]9(y-); (11) 



and the first power correction to transverse structure func- 
tion mm, 



f£\x b ,Q 2 ) 



~4ir 2 a s 




[ Q 2 


V 2iV c ). 



1 



< :, i B -^—T q {xB,A) (12) 



with correlation function, T q (xB, A), given in Eq. JSJ. 

The four-parton correlation functions T 9 's are nonper- 
turbative and must be taken from experiments. To es- 
timate the magnitude of T g 's, we could choose a simple 
ansatz [H?] 

T i {x,A) = X 2 A x l 5 <l> i/pA {x,A) (13) 

for i = q, g in terms of the corresponding twist- 2 effective 
nuclear parton distribution fa/A- We choose this form be- 
cause we expect the x-dependence of the probability to 
detect the hard parton to be essentially unaffected by the 
presence or absence of an additional soft scattering. In 
Eq. I |l 3(1 A is assumed to be a constant with dimensions of 
mass. This ansatz facilitates the comparison to data 0. 

If we further assume that a nucleus is made of a group 
of loosely bound color singlet nucleons, packed in a hard 
sphere of radius RA 1 / 3 , we can approximate the matrix 
element of nuclear state in Eq. |J5J into a product of matrix 
elements of nucleon states 1511201 . 



+ r 



{PA\m\ F 2 (0) qUDll-V 



(14) 



.4 
2^ 



IRA 1 ' 3 



(F+«F a +) (p\q(0)±- q (y-)\p) 
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(a) 




Fig. 5. Sketch for the ^-dependence of normal (a) and satu- 
rated (b) parton distributions. 



with p = pa/ A and the two-gluon matrix element |3] 
(F +a F a +) = i- / ^-(p\F+ a (0)F a +(y-)\p)0(y-) 



p + .1 2tt 



Substituting Eq. 1)15(1 into the definition of quark-gluon 
correlation function, T q (x, A) in Eq. J5J, we drive 



and 



37ra s 
8Q 2 i? 2 



<^ 4 



(16) 



/3 _ 



q A (x B ,Q Z ) (17) 



where g_A defined in Eq. JBJ| with p = pa/ A. In deriv- 
ing Eq. l|T7|. we used A 1 / 3, - 1 instead of A 1 / 3, for the 
line integral between the two pairs of field operators in 
Eq. I|15|l so that the nuclear effect vanishes for A = 1 [3]. 
From Eq. I(17JI and known generic x-dependence of parton 
distributions, as sketched in Fig. |S{a), the leading coher- 
ent power correction suppresses the DIS cross section or 
structure functions at small xb- 

Hxb is small enough, the virtual photon could interact 
coherently with all nucleons inside a large nucleus. From 
the Feynman diagrams in Fig. QJc) with 2N gluons, we 
derive the leading Nth coherent power corrections to the 
transverse structure function [3|, 



t T (x B , 



(F +a F a + ) (A 1/3 - 1 



3ira s 
8Q 2 R 2 

1 d N 

x 2£ e ^^^^(* B ,Q 2 ).(18) 

q B 



1 



Summing over all leading A^-enhanced power correc- 
tions, the first column of the right-hand-side of Eq. Q), 
we obtain (3] 



' A1LT W<3 2 



X xi 



d" 



dx% T 



(19) 




0.0001 0.001 0.01 0.1 0.001 0.01 0.1 



Fig. 6. From 0, all-twist resummed F2(A)/F2(D) calcula- 
tion versus DIS data on nuclei |25| . Data-Theory Ad-t is also 
shown. 



with a characteristic scale of the power corrections 

e = ^ (F+«F+ a ) . (20) 



gR2 \ ocf 

If we approximate N ~ oo in Eq. (|19|l . we obtain, 

F A (x B ,Q 2 ) « F^ (LT) (x B (l + Z\),Q 2 ) (21) 

with Zi ee (A 1 / 3 — 1) S, 2 /Q 2 , a shift in i^. Similar expres- 
sion was derived for the longitudinal structure function 

With only one unknown matrix element, (F +a F a + ), 
our calculated results can be easily tested and challenged. 
From Eq. (|21|) . the nuclear dependence in structure func- 
tions should come from two distinctive sources: universal 
A-dependence in leading twist parton distributions and 
process sensitive ^4-dependence from power corrections 2 
13] • By comparing our numerical results, evaluated with 
CTEQ6 PDFs EH, with the data, we extract the maxi- 
mum size of power corrections. For £ 2 = 0.09 — 0.12 GeV 2 , 
our calculated reduction in nuclear structure functions is 
consistent with the xb~, Q 2 - and ^4-dependence of the data 
as demonstrated Fig. El 

The predictive power of factorization approach resides 
in the universality of unknown matrix elements. Without 
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Fig. 7. Sample DIS diagrams have initial-state coUinear diver- 
gence involving parton recombinations. 



additional unknown matrix elements, our predictions for 
neutrino DIS cross section on an iron target are consistent 
with NuTeV data [22ES1- At small and moderate Q 2 , our 
prediction gives a very good description of SLAC E143 
data on B = <jl/&t of 12 C target |31I24"| 

The resummed power correction has a nontrivial de- 
pendence on the target A and kinematic variables xg and 
Q 2 because its connection to the ^-dependence of the ac- 
tive parton distribution. The effect disappears if the par- 
ton distribution is saturated at a very small sketched 
in Fig.OJb). The leading power correction calculated here 
is a consequence of the limited phase space in a given col- 
lision. If we would integrate over xb from — oo (or zero) 
to co for a purely inclusive process (or a process with an 

infinite collision energy) , the F^ N ^ vanishes for N ^ due 
to the unitarity condition [8J. 

QCD factorization systematically moves all coUinear 
divergences from partonic scattering into long-distance ma- 
trix elements. At the leading power, DGLAP evolved par- 
ton distributions are sufficient to absorb all collinear diver- 
gences in 7*-parton scattering cross section to ensure the 
infra-safety of all coefficient functions in the first row of 
Eq. 0). However, DGLAP evolved parton distributions do 
not remove collinear divergences involving multiple parton 
recombination, like those in Fig.|7Jb) and (c). A modified 
evolution equation taking into account the parton recom- 
bination slows down the fast growth of parton density at 
small x 20,20], and keeps a positive density of small- x 
gluons at low Q 2 [2Z|. All order resummed power correc- 
tions - coherent parton recombinations to parton evolu- 
tion should be vary valuable for approaching the region 
of parton saturation from the perturbative side and for 
getting a reliable picture of small- x partons |28U29j . 



4 Power correction to Drell-Yan cross section 

It was shown by the NA38 and NA50 Collaborations that 
muon pair production for dimuon invariant mass between 
4> and 3/tp, known as the intermediate mass region (IMR), 
in heavy nucleus- nucleus collisions exceeds the expectation 
based on a linear extrapolation of the p — A sources with 
the product of the mass numbers of the projectile and tar- 
get nuclei 130] . The excess increases with the number of 
participant nucleons, and the ratio between the observed 
dimuon yield and the expected sources reaches a factor 
of 2 for central Pb-Pb collisions. There have been a lot 
of effort to attribute such an excess to the enhancement 
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Fig. 8. Sketch for Drell-Yan process in hadron-nucleus collision 
with multiple interactions internal to nucleus (a) and sensitive 
to hard scattering (b); and the nuclear modification factor, 
Rab(Q), as a function of dilepton mass Q. 



of open charm production |31| . thermal dimuons produc- 
tion |32| . and secondary meson- meson scattering in nu- 
clear medium [33]. As shown in Ref. [30], the Drell-Yan 
continuum is the dominant source of the dilepton produc- 
tion in this mass range. An enhancement in the Drell-Yan 
continuum is much more effective than all other sources 
for interpreting the observed excess. 

In hadron-nucleus and nucleus- nucleus collisions, more 
partons are available at a given impact parameter. Before 
the hard collision of producing the lepton-pair, partons 
from different nucleons can either interact between them- 
selves, as sketched in Fig. [Sfa) which leads to the univer- 
sal nuclear dependence in parton distributions, or interact 
with the incoming parton, as sketched in Fig. [BJb) which 
gives the medium-size enhanced power corrections |34| . 
For the kinematics of IMR Drell-Yan process, parton dis- 
tributions have a week nuclear dependence and produce a 
small reduction to the cross section. On the other hand, 
the medium size enhanced power corrections are process 
dependent, and actually increase the production rate and 
become more important when Q 2 decreases [34j . Let in- 
clusive Drell-Yan cross section in nuclear collisions be ap- 
proximated as 



da AB 

dQ 2 



AB NN 



da 



(D) 
AB 



dQ 2 



dQ 2 



AB- 



. da 



(S) 
NN 



dQ 2 



1 + Rab(Q)] 



(22) 



where superscripts (S) and (D) represent the single and 
double scattering, respectively. The Rab(Q) defines a nu- 
clear modification factor of Drell-Yan continuum. Fig.[S^c) 
shows the calculated Rab(Q) from medium-size enhanced 
double scattering [33]. In evaluating Fig. [81c), we used the 
same quark-gluon correlation function, T q (x,A), and ne- 
glected the A-dependence in parton distributions. There- 
fore, we expect that the true enhancement to the Drell- 
Yan continuum might be slightly smaller than Rab{Q) in 
Fig. [Sfc), which is a very significant effect. 

It is perhaps surprising that the medium-size enhanced 
power correction to DIS and Drell-Yan cross sections car- 
ries a different sign. The sign difference is a consequence 
of the kinematic nature of the calculated leading power 
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(a) (b) (c) 

Fig. 9. DIS with a space-like hard scale 






t'p'-Ax'p' 

(a) (b) (c) 

Fig. 10. Drell-Yan with a time-like hard scale 



corrections. The sketch in Fig. [5] shows that the sum of 
multiple soft rescattering to the final-state parton in DIS 
effectively broadens the parton's transverse momentum 
(or gives the parton an effective mass [22]); and the 
on-shell condition of the final-state parton forces an ex- 
tra momentum fraction Axp to the incoming parton. On 
the other hand, the multiple initial-state soft rescattering 
in Drell-Yan process broadens the incoming parton's mo- 
mentum, and the kinematics (or the measured q^) of the 
observed dilepton pair forces the active parton from the 
target to take away some of its momentum, a negative 
shift in momentum fraction, as shown in Fig. HOf cl . 



5 Transverse momentum Qt broadening 

Coherent multiple rescatterings not only modify the pro- 
duction rate of inclusive cross sections, but also affect the 
momentum distribution of produced particles. Although 
the amount of broadening due to each soft rescattering is 
too weak a scale to warrant a reliable calculation, an av- 
eraged broadening in a hard collision could be a physical 
quantity calculable in pQCD |2). 

In Ref. |35| , Drell-Yan transverse momentum broaden- 
ing was calculated in pQCD by evaluating the lowest order 
soft rescattering diagram in Fig. Illf a). Although the di- 
rect Qt modification from soft rescattering to da/dQ 2 dQ T 
is not perturbatively calculable because of the size of small 
Qt kick, the Qt broadening, 



(Qt) 



dQrp Qrp 



da 



dQ 2 dQ\ 



dQ T 



2 lno\ 

dQ 2 dQ\ ' ( ' 



is calculable [2BH1- It was found that the Drell-Yan trans- 
verse momentum broadening in hadron-nucleus collisions 
can be expressed in terms of the same quark-gluon corre- 
lation function T q (x,A) (not its derivative) |35| . 

47r 2 a s \ E 9 e 2 q J dx'f m {x') T q (r/x\ A)/x' 



~ 1 ' ' E q e\ J dx'fg /h (x') f q (r/x>, A)/x> 

(24) 

where runs over all quark and antiquark flavors, e q is 
the quark fractional charge, and r = Q 2 / s, in terms of the 





GO 



Fig. 11. (a) Lowest order diagram that contributes to the 
Drell-Yan transverse momentum broadening in hadron-nucleus 
collisions; (b) Sketch for the lowest order contributions to the 
averaged di-jet momentum imbalance in photoproduction. 



lepton-pair invariant mass Q and hadron-hadron center of 
mass energy, y/s. 

Adopting the model in Eq. QI3[) . the lowest order Drell- 
Yan transverse momentum broadening in Eq. (|24|l can be 
simplified as 



(Q 2 



T/4/3 



47r 2 a, 



X 2 A 1 '\ 



(25) 



By comparing Eq. (|25|l to data from Fermilab E772 and 
CERN NA10 experiments j23S3, it was found that 
~ 0.01 GeV 2 , which corresponds to ££ Y ~ 0.04 GeV 2 
that is about a factor of 2 smaller from what was extracted 
from inclusive DIS data, where the leading twist nuclear 
dependence was not included. The difference might be 
caused by the difference between initial-state and final- 
state rescattering effects [2]- 

Momentum imbalance between two final-state jets in 
photoproduction off a nuclear target was calculated by 
evaluating the diagrams like one in Fig. Illf b). and was 
compared with Fermilab E683 data ^Hj) with an assump- 
tion that the momentum imbalance between two jets is 
approximately the same as the momentum imbalance be- 
tween two final-state partons. Again, the calculated mo- 
mentum imbalance was expressed in terms of the same 
quark-gluon and gluon-gluon correlation function |I9j . By 
comparing with the momentum imbalance data for jet 
transverse momentum pt > 4 GeV |38|. it was found 
[TU| that A^ ijct = 0.05 - 0.1 GeV 2 , which corresponds to 
£ 2 ~ 0.2 GeV 2 that is about a factor 2 bigger than what 
was found in DIS. As pointed out in Ref. |2], high order 
corrections to initial-state and final-state soft rescattering 
could be very different and significantly reduce the differ- 
ence of the nonperturbative parameter. 



6 Nuclear suppression in hadron-nucleus 
collisions 

Dynamical nuclear effects induced by soft rescattering can 
be studied through the ratio of particle production rates 
in hadron-nucleus and hadron-hadron collisions. However, 
there is no observed variable, like xb in inclusive DIS, 
that can directly measure the coherence length of the 
hard scattering. All incoming parton momentum fractions 
are convoluted over, and strength of the collision is mea- 
sured by Lorentz invariants, like Mandelstam variables, s, 



Jianwei Qiu: Rescattering effects in hadron- 



7 



t, and u. Unlike in DIS and in Drell-Yan, there are both 
initial-state and final-state soft rescattering in hadronic 
collisions, which could in principle lead to medium size 
enhanced power corrections. Since the hard scales, s, t, 
and u, in hadronic collisions are often much larger than a 
couple of GeV 2 , effect of the medium size enhanced power 
corrections to hard probes in hardonic collisions is in gen- 
eral less significant. However, in the most forward (back- 
ward) region, the invariant t (it) could be much smaller 
than the other invariants, so that the power corrections in 
1/t (1/u) could become very important |14j . 

Consider, for example, the single hadron inclusive pro- 
duction in hadron- nucleus collision as shown in Fig. ffi h). 
Once we fix the momentum fractions x a and z\ , the effec- 
tive interaction region is determined by the momentum 
exchange </ M = (x a P a — P c /zi) M . In the head-on frame of 
q — Pb, the scattered parton of momentum £ interacts co- 
herently with partons from different nucleons at the same 
impact parameter, just like that in DIS. Interactions that 
take place between the partons from the nucleus and the 
incoming parton of momentum x a P a and/or the outgoing 
parton of momentum P c j Z\ at a different impact parame- 
ter are much less coherent and actually dominated by the 
independent elastic scattering [HHI • Similar to the DIS case 
we find 51] that resumming the coherent scattering 
with multiple nucleons is equivalent to a shift of the mo- 
mentum fraction of the active parton from the nucleus in 

Fig. EDO, 



Xb 



x b [l + 



(-*) 



(26) 



with the hard scale t = q 2 = (x a P a — P c j z\) 2 and the color 
factor Cd depending on the flavor of parton d. C q ^ = 
1 and C g — Ca/Cf — 9/4 for quark (antiquark) and 
gluon, respectively. The shift in Eq. (|26J) leads to a net 
suppression of the cross sections, and the i-dependence of 
this shift indicates that the attenuation increases in the 
forward rapidity region. 

The top panels of Fig.ll2lshow the rapidity and trans- 
verse momentum dependence of R dAu (b) and R dA £ (b), the 
ratio of single and double hadron production, respectively, 
in minimum bias d+Au collisions. The amplification of the 
suppression effect at forward y\ comes from the steepening 
of the parton distribution functions at small Xb and the de- 
crease of the Mandelstam variable (— t). At high pti,PT2 
the attenuation is found to disappear in accord with the 
QCD factorization theorems 0Oj- The bottom panels of 
Fig. 1 121 show the growth of the nuclear attenuation effect 
with centrality. 

Dihadron correlations C2(A<p) = — associ- 
ated with 2^2 partonic hard scattering processes, af- 
ter subtracting the bulk many-body collision background, 
can be approximated by near-side and away-side Gaus- 
sians. The acoplanarity, Aip ^ tt, arises from high order 
QCD corrections and in the presence of nuclear matter - 
transverse momentum diffusion 39 . If the strength of the 
away-side correlation function in elementary N+N colli- 
sions is normalized to unity, dynamical quark and gluon 



i 1 i 1 

1.5 GeV, Min.bias 




Fig. 12. From |14| . Suppression of the single and double in- 
clusive hadron production rates as a function of pr in d + Au 
collisions at *J~s = 200GeV 2 at different rapidities (Top) and 
at different impact parameters (Bottom). The trigger hadron 



' • P + P \ STAR ' 

d+Au i f data 




A<t> [rad] 



A<(> [rad] 



Fig. 13. From |14) . Centrality dependence of Ci^Atp) at vari- 
ous rapidities and moderate (left) and small (right) transverse 
momenta. Central d + Au and p + p data from STAR |43| . 



shadowing in p + A reactions will be manifest in the at- 
tenuation of the area A Far = Rp A h2 {b) |T3|. 

The left panels of Fig. 1131 show that for p^ x = 4 GeV, 
Pt 2 = 2 GeV the dominant effect in C<2,{A(p) is a small in- 
crease of the broadening with centrality, compatible with 
the PHENIXg^l and STARg3] measurements. Even at 
forward rapidity, such as j/i = 2, the effect of power cor- 
rections in this transverse momentum range is not very 
significant. At small pr x =1.5 GeV, pr 2 = 1 GeV, shown 
in the right hand side of Fig. ^3 the apparent width of 
the away-side Ci^Akp) is larger. In going from midrapidity, 
yi=2/2 = 0, to forward rapidity, yi = 4, yi = 0, we find a 
significant reduction by a factor of 3 - 4 in the strength of 
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dihadron correlations. Preliminary STAR results |44| are 
consistent with our predictions. 



7 Summary and outlook 

We have briefly reviewed a pQCD factorization approach 
to coherent multiple scattering and argued that the A 1 / 3 - 
type medium size enhanced power corrections caused by 
multiple soft rescattering can be consistently calculated 
in pQCD. We presented calculations of rescattering ef- 
fects in inclusive cross sections as well as transverse mo- 
mentum broadening (or the moment of particle transverse 
momentum distributions.). By studying the coherent mul- 
tiple scattering, we can probe new sets of fundamental 
and universal multiparton correlation functions in nuclear 
medium. These new functions provide new insights into 
the nonperturbative regime of QCD. 

Initial pQCD calculations of multiple rescattering ef- 
fects is very successful in understanding nuclear depen- 
dence of inclusive cross sections. For a wide range of ob- 
servables, the estimate for the only unknown parameter is 
reasonably consistent. With a single unknown parameter, 
the calculations describe many observables and their nu- 
clear dependences fairly well. The factorization approach, 
with its intuitive and transparent results, can be easily ap- 
plied to study the nuclear modification of other physical 
observables in p + A reactions. The systematic incorpora- 
tion of coherent power corrections provides a novel tool 
to address the most interesting transition region between 
"hard" and "soft" physics in hadron-nucleus collisions. 

This work is supported in part by the U.S. Department 
of Energy under Grant No. DE-FG02-87ER40371. 
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